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Abstract 

The propagator for gluons in a space-time dependent field is derived. This is accomplished 
by solving the equation of motion for the gluonic Green's functions. Subsequently a relationship 
between the quark and the gluon propagator is presented. With its help characteristics of the 
bosonic propagator can be obtained from those of the fermionic and vice versa. Finally, this 
relation is discussed for the special case of ultrarelativistic collisions in the semiclassical limit. 
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I. INTRODUCTION 



In Refs. 0, 0, HI the full retarded fermion propagators in the presence of arbitrarily space- 
time-dependent classical fields were studied. Here the generalisation to the corresponding 
gluon propagator is presented. Correlators of this kind are, for example,needed in order to 
study the production of the respective particles by vacuum polarisation 0, 0, IE 0, IE EL IE HI 
or calculate induced currents and condensates. Such phenomena are of importance for the 
physics of the early universe |9| and ultrarelativistic heavy-ion collisions together with the 
quark-gluon plasma (QGP) [10(. The strong classical fields are the common feature in these 
systems, while the detailed characteristics of the field configurations differ. After provid- 
ing the general formalism, this paper focusses on particle propagation in ultrarelativistic 
collisions. 

A lot of effort is made to study the QGP's production and equilibration in nuclear 
collision experiments at the Relativistic Heavy-Ion Collider (RHIC) and the Large Hadron 
Collider (LHC), currently under construction. It is a widely used assumption that the ini- 
tial state in heavy-ion collisions is dominated by gluons. Due to their large occupation 
number and as a first approximation one can begin the system's description with a classi- 
cal background field. Quantum fluctuations have to be investigated subsequently. Those 
considerations have also lead to the McLerran-Venugopalan (MV) model 2|- There, the 



projectiles are represented by classically interacting colour-charge distributions on the two 
branches of the light-cone. One observes that quantum corrections are enhanced by kine- 
matical factors. They can be included into the distributions by means of the JIMWLK 
equation 1^ . 



In any case, the high occupation-number bosonic fields are that strong that multiple 
couplings to the classical field are not suppressed even for a small coupling constant g. In 
the absence of other scales they have to be taken into account to all orders. The lead- 
ing quantum processes concern terms in the classical action of second order in the fields 
(fermions, antifermions, and bosonic quanta). From those terms, the inverse of the respec- 
tive two-point Green's functions can be read off. Their solutions give the propagators of 
the corresponding particles to all orders in the classical field. Depending on the boundary 
conditions that are imposed, the correlator will serve to obtain expectation values (in-in- 
formalism, e.g., retarded and advanced propagator) or probabilities (in-out-formalism, e.g., 
Feynman propagator). In the following, the first possibility will be pursued mostly. 

The paper is organised as follows: Section II contains the derivation of the equation of 
motion for the gluon Green's functions within the framework of the background-field method 
of QCD. Section III derives the retarded gluon-propagator as a solution of the equation of 
motion in background-field Feynman- gauge. Section IV gives the retarded Green's function 
solution for the quadratic Dirac-operator. In section V the interrelation between the retarded 
gluon propagator and the retarded quadratic Dirac-propagator is worked out. Additionally 
the connection to the retarded linear Dirac-propagator is shown. These results are used 
in order to translate the known characteristics of the two-point function for the first-order 
Dirac-equation into those of the other two correlators. Section VI applies the findings of the 
previous sections to ultrarelativistic heavy-ion collisions. The cases of nucleus-nucleus and 
nucleus-nucleon collisions are addressed as well as that of the propagation of gluons with a 
large momentum scale. In the last situation all terms for the retarded propagator are spelled 
out and the Feynman propagator is available, too. Section VII summarises the results. 

Throughout the paper the metric tensor is given by: g^ u = diag(l, —1, —1, —1), angular 



2 



momenta are measured in units of h, and velocities in fractions of the speed of light c. From 
hereon, the coupling constant is included in the definition of the classical field: gA^ ld = A^ ew . 
v represents the three- vector of the spatial components of any four- vector v. The convention 
for Fourier transformations of one-point functions is: 



f(x,y)= I ,i\ ,t\ e->°-'e +i ™f(q,p). (2) 



that for two-point functions: 

?A = ( 

(2tt) 4 (2tt) 4 

II. THE EQUATION OF MOTION FOR GLUON GREEN'S FUNCTIONS 

In order to derive the requested equation of motion for the gluonic Green's functions, let 
us start with the classical action for the bosonic sector of the system: 



1 
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with the Lagrangean density for the gauge field 



S = \J d 4 xC G (3) 



C G = ^F a :F a (4) 

where the field tensor is defined as, 

F lu ,=i[D li ,D v ] (5) 

with the covariant derivatives: 

Dfi = <9 M — (6) 

The field tensor is explicitely given by: 

F; u [A] = d,Al - d v A; + r hc A\Al. (7) 

The splitting of the gauge field A^ into a classical A^ and a quantum field: 

Afj, = A^ + Q„, (8) 

allows to express the field tensor as a functional of two fields: 

F; V [A + Q] = F; v [A\ + Df[A]Q b u - Dt b [A]Ql + f abc QlQl- (9) 

-D^ 6 [A] stands for the covariant derivative as a functional of the classical field A^ only. 
The Lagrangean density now is composed of a sum of terms containing different powers of 
the quantum fields. The terms of second order determines the equation of motion for the 
propagator of the quantum fluctuations in the classical background. It is given by: 

4? = \ {{Df[A\Ql){D ac »[A\Q c ») - (D; b [A]Q b u )(D a ^[A]Q c n - f abc F« u Q b »Q c »} . (10) 
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Additional terms of second order in the quantum fluctuations arise from the gauge-fixing 
term. The gauge fixing term to be added to the Lagrangean density reads for the background- 
field Feynman- gauge: 

C GF = L% = - l -{Df[A}Q b »){Dl c [A]Q™) (11) 

As the Lagrangean density is always integrated over in order to obtain the action, the 
expressions for the second order terms can be transformed by virtue of partial integrations. 
The representation as Lagrangean density is kept for the sake of simplicity: 

Cf = l -Q b » {D b x a [A]D acX [A] 9flu - Dfr[A]D?[A]} Q cv (12) 



Cf p = ±Qf«'D*[A]D?[A](r. (13) 

Finally, the differential operator for the equation of motion for the gluonic Green's func- 
tions can be read off by omitting the quantum fields. Thus the equation of motion reads: 

T-\x)T, u (x,y)=6 (4 \x-y)g tlu (14) 

with the differential operator (l8| : 

r _1 = D[A) ■ D[A). (15) 

Further the contribution to the Lagrangean density originating from the Faddeev-Popov 
ghosts in background-field Lorenz-gauge is given by: 

C FP = -{Df[A\ X *){D™[A + Q]x c ) (16) 

Partial integration yields: 

C FP = x b ^D b ;[A}D a ;[A + Q} X c (17) 
Therefore, the equation of motion for the ghost propagator in the classical field reads: 

T- 1 (x)T(x,y)=5^(x-y). (18) 
Note the connection between the quantum gluon and the ghost propagator in this gauge: 

IV = Tg^, (19) 
whence they can be used as synonymous to each other for most practical purposes. 
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III. THE RETARDED GLUON PROPAGATOR 



First look at the homogeneous equation for the gluon (ghost) propagator: 

{D[A(x)]-D[A(x)}}T H (x,y)=0 (20) 

Explicitly one has: 

{d 2 — 2iA d — A + 2iA ■ d — i(d ■ A) — A - A}(x)T H (x,y) = (21) 

where the time derivatives have been singled out, because the boundary conditions for the 
retarded propagator are given for equal time. The differential operator is taken at the space- 
time point x, which is only denoted once for the sake of brevity. A stands for the Laplace 
operator. In the term (d ■ A) the derivative does not act accross the rounded brackets. The 
previous equation can be reexpressed as 

f-\x)f H (x,y) = (22) 

with 

f- 1 = a d + f 1 (23) 

where 

f" 1 = —iuoAo — a + + cr_{— A + 2iA ■ d — i(d ■ A) - A ■ A} + ia 3 A (24) 

while 

2f H = a (r« + d T {2) ) + 2a + rg } + 2^d T^ + a 3 (T { ^ - d T%) (25) 

where r# and are two independent solutions of Eq. (j2*UJ) . Oj with j G {1, 2, 3} are the 
Pauli matrices and a the corresponding unit matrix. a± are defined as 2cr ± = cr 1 ± i<j 2 - Eq- 
(1221) is solved by: 

t H (x,y) =Pexp|£°^of- 1 (eo,x)}. (26) 

Choosing a different boundary condition than Tn(x,y) = cr at xq = yo would only lead 

to overall factors to the right of the path-ordered exponential, which, for a homogeneous 

differential equation does not lead to independent solutions. 

The retarded gluonic propagator is to vanish for negative time-differences: T(x, y) = for 
~ Ho < 0. Further, it is to be continuous at x = y and there it is to have a discontinuous 

first temporal derivative: \im xo ^ yo+0 d (x)T(x,y) — lim xo ^ yo _ d (x)T(x,y) = 5^ 3 \x — y). 

Hence, the retarded propagator can be expressed as: 

T R (x,y) =f {2 \x,y)6(xo-y )6 i3 \x-y) (27) 

Like in 0, |j| this solution is fit to be analysed in various configurations with the help of the 
general resummation formula (j3*7jl and the group property ()38|) . However, here information 
on the characteristics of the gluonic two-point function are to be obtained by comparing it 
to quark correlators. 
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IV. THE RETARDED QUADRATIC FERMION-PROPAGATOR 

Starting out from the linear Dirac-equation 

{z 7 • D[a(x)] - m}G(x, y) = 5 {4) (x - y) (28) 

in the presence of the gauge field a in fundamental representation, the quadratic differential- 
equation is obtained through the substitution 

G(x, y) = {-ry ■ D[a(x)\ - m}g(x, y) (29) 

and one finds 

{D[a} -D[a] - a:F/2 + m 2 }(x)g(x,y) = 5 {4 \x - y) (30) 

where use has been made of the definition of the field tensor (jSJ) and with the commutator 
of the 7 matrices 2a fJ ' u = — i[7 M , 7 1 ']. 

In analogy to the gluonic case, the homogeneous part of this equation is equivalent to: 

g-\x)g H (x,y) = 0, (31) 

with 

g- X = vofo + r X (32) 

where 

g^ 1 = —ia Q a — a + + cr_{— A + 2ia ■ d — i(d ■ a) — a ■ a — cr:F[a}/2 + m 2 } + itT 3 a (33) 
while 

Tg H = M9h + dogP) + Zv+gP + ^d g^ + a 3 {g$ - d gP) (34) 

where g$ and g^ are two independent homogeneous solutions of Eq. (J30J). Then Eq. 
is solved by: 

9h(x,v) ="Pexp|y d^ g~\^o,sc)j (35) 

The boundary conditions in the fermionic case are the same as in the 
bosonic, i.e., g(x,y) = for x - y < 0, lim xo _> yo+0 g(x,y) = lim xo ^ yo _ g(x,y), and 
lim xo ->y 0+0 d (x)g(x,y) - linxro-^o-o d (x)g(x, y) = 5 (3) (f -y), whence the retarded propa- 
gator can be expressed as: 

9r(x,v) = g i H ) (x,y)6(x - y )6 {3) (x- y). (36) 
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V. INTERRELATION BETWEEN THE PROPAGATORS 



With the help of the general resummation formula derived in j^: 

Pexp{£Vo[£(fo) + <?(&)]} = 

= Pexp{jfVo£(fo)} x 

r%o ( ryo •) ( r£o 

xPexp / d£ Vexp\ dz B(z )\ C(£o)P exp I dz B(z ) 
\Jyo > {Jyo 

= Pexp / (i^exp / dz B(z )>C(^o)'Pexpi dz B(z ) 

xPexp{£Vo5(£ )}, (37) 



which is based on the group property valid for path-ordered exponentials: 

7>exp{£°d£o5(£o)} = Pexp {jT d&flfo)} x Pexp |£° d£ £(£o)} , (38) 

the bosonic propagator can be identified in the fermionic. This is achieved by making the 
choice C = c = a_{—cr:F[a\/2 + m 2 }. The homogeneous solution of the quadratic Dirac- 
equation can then be reexpressed as: 

9H(x ,yo;x) = Pexp|+ J d^ f' H (x ,^x)c^o,x)t' H ^ ,x ]x)^t' H (x ,y ]x). (39) 

The y in all the arguments has been left out, because due to the 8^(x — y) present in the 
propagator, it is always taken equal to x in the end. The bosonic function has to be taken 
in the fundamental representation. This supersymmetric replacement has been indicated by 
the prime. 

Ultimately, the inverse relation is required, whence the task of identifying the fermionic 
propagator in the bosonic arises. To this end, start out with the bosonic homogeneous 
solution (|26|) . First multiply with unity in the space of the 7 matrices everywhere. The 
inverse operation is taking one quarter of the trace over the elements of the Clifford algebra, 
which has not been denoted seperately for the sake of brevity. Afterwards, the application 
of the resummation formula (}3*Tj) with C = — c results in: 

^H(x ,y ;x) =Pexp j-^ d^g' H {x G , £ ; x)c'(£ , x)g' H (£ , x ; x)| g' H {x , y ; x) (40) 

The primes indicate that the fermionic functions are constructed in the adjoint representa- 
tion. The inserted matrix d can be split into the contributions from the mass and the field 
tensor, respectively. Resumming all terms involving the mass corresponds to replacing every 
massive quark function g' H by a massless one g' h : 

t H (x ,y ;x) ="Pexp d£ g' h (x , £ ; x)cr:F[A(£ , x)]cr_g' h (£ , x ; x)J g' h (x , y ; x) (41) 



7 



Again this expression is fit to be studied with the help of Eqs. (|3*7j) and (|3H|) along the 
lines of Refs. 0, |j| . Here we shall concentrate on an expansion in the insertions of the field 
tensor. To lowest order this yields for the relevant cr + component: 

f S 1(2) (x , y ; x) = gf ] (x , y ; x) . (42) 

The corrections arising in higher orders of the expansion include further couplings to the 
field tensor. The a + component of the first order reads: 

F l H {2 \xa,yo;x) = \ I d£ gP(xo,£o;x)a:F[A(£o,x)}gf\£o,y ;x), (43) 
* Jyo 

where use has been made of the relations u±u T cr± = <r± and (cr ± ) 2 = 0. In fact, to all orders 
only the a + component of g' h contributes to V K H ' . Multiplying the previous two equations 
with 9(x — yo)S^(x — y) in order to obtain the contributions to the retarded propagator 
(1271) leads to: 

T [ S(x,y)=g' r (x,y) (44) 

and 

r£W) = \J d 4 U(x,0^[A(0}g'M,y), (45) 

where 9(xo — yo) = 9(xo — £o)#(£o — ^o) for £ [xq, yo] and afterwards the definition of the 
massless retarded quadratic fermion-propagator g r in accordance with Eq. (j^fijl have been 
exploited. In vacuum configurations of the gauge field Eq. (|44jl becomes exact and the two 
propagators are identical. Last but not least, the retarded quadratic fermion-propagator gu 
can be constructed from the retarded propagator Gr for the linear Dirac-equation. Here the 
massless case is needed: 

g r (x, y) = —J d 4 zG r (x, z)G r (z, y). (46) 

which provides the link to the results presented in Refs. (HQ, El- The minus sign originates 
from the sign of the covariant derivative in Eq. (|29)1. 



VI. SEMICLASSICAL ULTRARELATIVISTIC HEAVY-ION COLLISIONS 

The semiclassical description of an ultrarelativistic heavy-ion collision is usually begun 
with a current of colour charges moving along the light-cone 0|. In order to obtain the 
classical gauge field the Yang-Mills equations have to be solved in the presence of this current. 
In Lorenz gauge the field consists of two different contributions. One of them is the radiation 
field inside the forward light-cone. The other is made up of the Weizsacker- Williams (WW) 
sheets. These correspond to the Coulomb fields of the different colour charges boosted 
into the closure of the Lorentz group. Therewith they become ^-distributions in the light- 
cone coordinates. This approximation is justified for the description of particles with low 
longitudinal momentum (mid-rapidity). On the one hand they have a low resolution in this 
direction. On the other they cannot be comovers of the charges on the light-cone. For 
comoving particles the plane- wave expansion would not be applicable 0. 



S 



In the Lorenz gauge for the classical field, the WW contribution to the gauge-field takes 
the form: 



.4 



WW, 



A 



WW, 



A 



WW, 



X) 



X) 



X) 



E *a(tfl)n £ <5 
n L =l 



N L 



g_ 

'2n 



Nr 



E UWnJ 



X- 



Xa 



n R =l 



(&5 



IriR 



In A 



In A 



x T - (6r 







(47) 



A is an arbitrary constant, regularising the logarithm, which does not appear in quantities 
like the field tensor. The t a are the generators of SU(3) C . The (t^ R ) n represent the colour 
of the charges. 

In general, the presence of a second nucleus leads to a precession of the charges of the first 
nucleus and vice versa which is a manifestation of the covariant conservation of the current. 
This causes further sheet-like contributions to the gauge field which can be combined with 
the WW fields by modifying the charges. Up to the next-to-leading order in perturbation 
theory they read fl3^ : 



ita)n L 



fta )"r 



Nr 

(t L a )t- a s E UM)TM)T B 

n R =l 
N L 

+ E fa b c{t L b )Zitc)nR d 

n L =l 



- (6j) nH hiA|x T - (&£) n J + 0{a s 



X- 



{b L _) nL \n\\x T -(b^) nL \ + 0(a s 2 ). 



(4f 



with the initial colours {ta ,R Yn L R - To all orders, the modifications amount to Wilson 
lines over the gauge field along the branches of the light cone. Hence, the preces- 
sion terms are absent in adequate, i.e., light-cone gauges. Colour neutrality requires 
YlnL=\ ta(ta)nL = = Y^nR=i^aita) % n R - Higher-order perturbative calculations for the gauge 
field as solution of the Yang-Mills equations are only of value in the presence of hard energy- 
or short time-scales. Up to now, non-perturbative solutions for this problem have been 
obtained only in transverse lattice calculations p3 |. In any case, the general form of the 
field — continuous inside the forward light-cone and singular but integrable on the forward 
light-cone — allows to express the retarded propagator by a finite number of addends. 

The final result for the massless retarded fermion propagator G r for the linear Dirac- 
equation in the ultrarelativistic limit reads [l|: 

zG r (x,y)7° = +G h (x,y)6^(x -y)9{x Q -y ), (49) 

with the homogeneous solution of the Dirac-equation: 

G h (x,y) = G r h ad (x,y) + 

+ J d^Gl ad (x,0[T L (0^-) +T R (tm + )]GH d (Z,y) + 

+ J rf 4 eAGr^,OT L (0^(e-)Gr(e^)T^(r ? )5(r /+ )Gr(^2/) + 
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where for an even number of scattering centers in each projectile the scattering matrices are 
given by 

T L (x) + p + = p + V I] ex P H<(x T )], 

n L =l 
Nr 

T R (x) + p- = p-V J] ex P H<(f T )], 

n R =l 



with 



2na 



L,R 



-gt a (t L a 



' - In A 



xt — (b T 



™L,R 



" n L,R »"»\"a J n L,R 

while the homogeneous solution in the radiation field is given by 

i yo 



Gl ad (x, y)=Vex V U d^^d^x) + 7 • x)} 



(51) 



(52) 



(53) 



Higher orders in the T cannot contribute to Eq. (jBHj) . because lines of constant a;_ or x + 
can only be crossed once [H S| • 

This remains also true for the convolution of two retarded propagators in Eq. (J46j) . There 
are at most contributions with one coupling to each scattering matrix in g r (x,y). 



As already mentioned above, in nucleus-nucleus collsions the radiation field A rad is only 
known numerically up to now j 141 ] . This problem can be circumvented in proton- nucleus 
collisions by regarding the charge density of the proton as small ^5[. The Yang-Mills 
equations can be solved to all orders in the charge density of the nucleus and to the first 
order in the charge density of the proton. |2(j The propagator in the radiation field is replaced 
by its perturbative expansion up to the first order in the thus obtained radiation field. The 
scattering matrices for said case are obtained from the above expressions by keeping the 
nuclear scattering matrix T L (x) as it is and by replacing the one for the proton by the 
lowest non-trivial order: 

Nr 

T R ] {x) = -ip~ £ a R R (x T ). (54) 

n R =l 

Note that the addend where the radiation field and the protonic scattering matrix contribute 
simultaneously, is discarded, because it is of second order in the charge density of the nucleon. 



In the case where a large momentum scale is assigned to the propagated particle, be it in 
nucleus- nucleus or in proton-nucleus collisions, the equations simplify further ja, [l3j . Also 
the other scattering matrix in Eq. (|51|) is reduced to its lowest order: 

N L 

T^{x) = -ip + <(*t). (55) 

n L =l 

The Yang-Mills equations are only solved perturbatively to the first order in each of the 
charge densities to yield the radiation field A™ d . The propagator in the radiation field 
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is expanded to the first order in the field. At the end, only couplings to the sheet-like 
contributions or to the radiation field are kept. Taking them into account simultaneously, 
leads to higher order terms. Thus the retarded linear Dirac-propagator to this order is given 
by Eq. (JHJ) with: 

Gff(s,y)= G%(s-y) + 

+ J d^G° h (x - EffigMg-j + JgCOfg+j - *7°7 • <(fl]G£tf - y) + 

+ 1 AA^x - - i7)lgfa)*fo + )G£fo - y) + 

+ / d^ V G° h (x - 0T«(05(Z + )G° h (t - ry)^ (v)Kv-)G° h (v - v)- (56) 

G\(x — y) stands for the free homogeneous solution of the linear Dirac-equation without 
mass. Putting this into Eq. (|4*B|) an d subsequently into Eq. (jUj) yields: 

g?(x,y) = -g° r (x-y)- J d 4 zG° r (x - z)G] l (z,y) - J d A zG?(x, z)G° r (z - y) + 
+ J d^d 4 V G° r (x - - r?)!]^)^)^ - y) + 

+ y d 4 eAG°(a; - " v)T^(v)5(v + )G° r (v - y) (57) 

with 

9° r (x,y) = -J d 4 zG° r (x-z)G° r (z-y) (58) 

For the use with Eq. ()45|) the addends which are not underlined will not contribute at 
all. In connection with the radiated field tensor merely the free propagator in Eq. (|5l)Jl 
is needed. Hence one finds, still with some higher order terms to be removed in the first 
addend: 

= \f*Z£{x-e)a:F[A ww (M{Z-y) + 
+\ J d 4 U(x, QaiFlAWitMit - y) + 
+\ J d^g° r (x - Oa:F[A ww (£M(t, y) + 
+\ J d A tg° r (x - t)<r.F[Atf(t)]g° r (t - y) (59) 

with: 

y,(x,y) = + J ^(i-OK(W+^(0«(e-j) + 

+ / d^G° r ( X - + Tftmz+Mit - v)- (so) 

To this accuracy, the last contribution comes from the second order term of Eq. (j41|) : 
r [ $(x,y) = \jd^d 4 V g Ax-0o-:F[Ay w (0]g^-v)°--F[A w ^ 

+ \ J d^d 4 V g° r (x - Oor:F{A ww (0}gr(Z ~ rj)a:F[A™ (v)}g r (r] - y), (61) 
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where only the field tensors as functionals of the sheet-like configurations (|47|) 1 



F +T [A 



ww-i 9 



f E t«(ta)n L 5 [X- - (b L _) nL 
271 n L =l L 



XT ~ Q>T)n L 



and 



F- T [A WW ] = Jj- g * a (t*) nii <5 [x + - (6?) nH 



E 



XT - (&T)n L 
XT ~ (&£) 



II R 



x T - 



(62) 



(63) 



and all other components equal to zero, contribute. 

Here, as at most two field insertions of any kind per addend are allowed, the corresponding 
Feynman propagator can be obtained by replacing the free retarded propagators by free 
Feynman propagators. 



Gauge independent results can be obtained from an operator in different ways.for ex- 
ample, by multiplying with adequate gauge links or by averaging over all gauges The 
MV-like models fall back on averaging over an ensemble of charge distributions with a gauge- 



invariant weig ht E3. This averaging can be carried out equivalently over the WW fields 

0. 

In the case of a gluon propagating at a large momentum scale, the radiated field is a 
functional linear in A^ w and A^ w . Thus the same holds for the gluon propagator in 
this situation. Due to the colour neutrality of each projectile, the average over the charge 
densities or equivalently the WW fields are to vanish (A?± w j = 0. Therefore, the average 
propagator to this order coincides with the free propagator. The deviations from the free 
case start but at the next order or if the average over an operator is taken that contains 
higer powers of the propagator. 

Along the same lines, the propagator in the proton-nucleus case is, per definition of the 
approximation, a functional linear in the charge density or the WW field of the proton. This 
does not change if terms beyond Eq. (|4"5Jl are taken into account. Hence, after taking the 
average over the proton imposing colour neutrality, only the average over the propagator in 
the field of the nucleus is left. As before, deviations originating from the presence of the 
nucleon only start playing a role in the following orders or in averages over powers of the 
propagator. 



VII. SUMMARY 

Solutions for the retarded gluonic propagator in background-field Feynman gauge and the 
retarded propagator for the quadratic Dirac-equation have been derived to all orders in the 
classical field. Imposing the corresponding boundary conditions has singled out the retarded 
correlators. Subsequently, relations expressing each of these propagators as a functional of 
the other have been determined. The one constructing the bosonic two-point function from 
the fermionic has been appended with another one, linking the quadratic fermion propagator 
with the linear one. With the help of said equations, the results for Green's functions of the 
linear Dirac-operator from Refs. 0, 0, II can be and have been transferred to those of the 
quadratic and the bosonic operator. 
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This fact has been exploited for the case of ultrarelativistic collisions. For nucleus-nucleus 
collisions all pieces of information, which do not require numerical calculations for the ra- 
diated field have been exposed. For nucleon-nucleus collisions, one can go further, because 
they allow for an expansion of the Yang-Mills equations in the charge density of the proton. 
The building blocks for the description of this system have been provided. 

As a third system, the propagation of a gluon with a large momentum scale is investigated. 
In this case the Yang-Mills equations can be expanded in the densities of both projectiles. 
All possible terms up to first order in any of the charge densities have been derived. In this 
case, also the Feynman propagator can be obtained. At this order, averaging the propagators 
over the charge distributions along the line of the MV model, gives the free propagator. 
Deviations arise when averaging over operators containing powers of the propagator at this 
order or over the propagator at an higher order. In the case of a proton-nucleus collision, 
the averaging procedure leads to the average over the propagator in the field of the nucleus. 

It will be interesting to apply the presented correlators to describe phenomena like vacuum 
polarisation and quantities like induced currents or condensates in the gluonic case. Such 
investigations will be presented elsewhere. 
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